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The consolidation of colloidal particles in drying colloidal dispersions is influenced by various 
factors such as particle size and shape, and inter-particle potential. The capillary pressure induced 
by the menisci, formed between the top layer of particles in the packed bed, compresses the bed of 
particles while the constraints imposed by the boundaries result in tensile stresses in the packing. 
Presence of flaws or defects in the bed determines its ultimate strength under such circumstances. 
In this study, we determine the asymptotic stress distribution around a flaw in a two dimensional 
colloidal packing saturated with liquid and compare the results with those obtained from the full 
numerical solution of the problem. Using the Griffith's criterion for equilibrium cracks, we relate 
the critical capillary pressure at equilibrium to the crack size and the mechanical properties of the 
packed bed. The analysis also gives the maximum allowable flaw size for obtaining a crack free 
packing. 

PACS numbers: 82.70.Dd 



I. INTRODUCTION 

Dried colloidal particle films find use in a number of 
applications such as tapes for photography and magnetic 
storage [1], porous coated printer papers, coating vitamin 
tablets, synthetic opals, photonic crystals [2] etc. The 
macroscopic properties of the film such as its thickness, 
particle packing and the mechanical strength arc influ- 
enced by the drying rate, interparticle potential, particle 
size and shape and the modulus of the particles. When 
a dispersion of colloidal particles is dried, the particles 
concentrate, eventually reaching a close packed concen- 
tration. The liquid menisci on the top layer of particles 
compresses the packing while the substrate resists trans- 
verse deformation of the packing. Consequently, trans- 
verse tensile stresses develop in the packing and when 
these stresses exceed a critical value, the packing cracks 
resulting in a variety of crack patterns. Such cracks oc- 
cur not only in thin films such as paints and coatings but 
also in thick systems and over geophysical length scales 
such as in the case of dried river beds. 

Most of the experimental investigations of the crack- 
ing phenomenon in drying colloidal dispersions have fo- 
cused on the thin film geometry where stresses have been 
measured using the classical cantilever bending technique 
[1, 3, 4]. These measurements show that thin films of 
monodisperse colloidal dispersions containing identical 
particles crack at a critical stress that is independent of 
the particle size but varies inversely with the film thick- 
ness [3, 5-7]. In almost all cases, the film nucleates multi- 
ple cracks with crack spacing that varies linearly with film 
thickness [3, 8]. Experiments also suggest that irrespec- 
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tive of particle size or moduli, each dispersion has a max- 
imum crack free thickness below which the films do not 
crack. The critical cracking thickness is found to increase 
with particle size and moduli in the case of hard polymer 
and metal oxide particles. A number of investigations 
have also focused on cracking in confined geometries such 
as capillary tubes where the dispersion dries from one 
end resulting in a compaction front of packed particles. 
While the studies in this geometry have mainly focused 
on crack tip velocity and its relation to the speed of the 
compaction front [9-11], it is only recently that Dufresne 
and co-workers [12] have been successful in imaging the 
stress variation near the tip of a propagating (interface) 
crack at the interface of an elastomer and saturated col- 
loidal bed and extract the stress intensity factor from it. 
The stress decays as inverse square root of the distance 
from the crack-tip and is in line with the predictions of 
classical linear theory for fracture in brittle materials. 

On the theoretical front, Routh and Russel [13] have 
derived a constitutive relation relating the macroscopic 
stress to macroscopic strain in a drying film. They con- 
sidered the viscoclastic deformation of a pair of identical 
particles due to contact and interfacial forces and related 
the strain at the particle level to these forces. Next, they 
volume averaged the forces over all orientations to arrive 
at the macroscopic stress versus strain relationship for a 
drying film. In the absence of particle-solvent interfacial 
tension, the expression for the macroscopic stress tensor 
[3] for identical elastic spheres reduces to 
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where, e^- is the macroscopic strain, P is the capillary 
pressure, 4> rcp is the random close packing concentra- 
tion, G is the shear modulus of the particles and M is 
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the number of contacting neighbors. The constitutive 
equation is an improvement over the traditional poroe- 
lasticity models [14-16] as the former accounts for the 
nonlinear deformation at the particle level and the influ- 
ence of particle size, modulus and packing characteristics 
on the macroscopic deformation field. The model has 
been successful in predicting not only the stress profile 
in drying films of both film forming and cracking systems 
[17], but also in predicting many aspects of the cracking 
mechanism in the latter [3]. More recently, Russel et. 
al. [18] have improved on the above relation by adopting 
the Hertzian contact mechanics at the particle pair level. 
The final constitutive relation is also non-linear with the 
stress varying as three halves power of the strain. Using 
this relation, they determine the capillary pressure nec- 
essary either to open an infinite crack in a flawless film 
or to extend pre-existing flaws of finite lengths. Their 
results suggest that flaws which are a fraction of the film 
thickness are sufficient to initiate cracks that would prop- 
agate across the sample at pressures modestly greater 
than obtained from the energy argument. In a related 
study Man and Russel [19] demonstrate experimentally 
the role of flaws in nucleating a crack and show that the 
critical stress obtained from the energy argument only 
gives the lower bound. 




FIG. 1. A crack of length 2a embedded in a packed bed. The 
bed is stressed in '1' direction. We consider the case where 
a < h. 



In this study, we determine the stress field near a crack 
tip along with the shape of the crack that is present in a 
two dimensional particle packing saturated with solvent. 
The flaw is embedded inside the colloidal packing and 
the size of the flaw is much smaller that any other dimen- 
sion of packing, say film thickness (Figure 1). Further, 
when the crack dimension in the out of plane direction is 
larger than a, then the only length dimension relevant to 
the problem is a, and the situation becomes amenable to 
plane stress/strain analysis [20]. A packed bed made of 
an array of colloidal particles can be considered to be a 
collection of polycrystalline aggregates with pre-existing 
flaws. These flaws may be attributed to micro-cracks, 
grain boundaries between the clusters of ordered pack- 
ing of mono-dispersed particles, dissimilar pores inside 
the colloidal bed etc. Nucleation of a crack under these 
circumstances changes the stress field close to the crack 
with stress concentration at the crack tip. In this work, 
the stress and strain fields are linearized about the pre- 
crack state to determine the disturbance displacement 
field immediately after the opening of a mode-I crack. 
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FIG. 2. A crack of length 2a embedded in a packed bed. 
The bed is stressed in '1' direction. Shaded region shows the 
region over which the analysis is performed. T{ is the traction 
along crack surface. 



These results also yield the stress intensity factor for 
the two dimensional elastic field which is then related 
to the surface energy using the well known Griffith's cri- 
terion for equilibrium cracks. The calculated quantities 
are then compared with the numerical solution for the 
full problem. The calculations show that the dimension- 
less critical capillary pressure required to open a crack 
varies inversely with the crack length to the two thirds' 
power and depends on a dimensionless parameter that 
measures the ratio of the elastic to surface energy. A 
simple scaling analysis reveals the essence of the results 
to follow. Since a ~ Ee° 2 , where '-E' is effective modu- 
lus of the packing and e° is the characteristic strain in 
the packing, the elastic energy recovered on the opening 
of a crack of length 'a' in a packing of unit thickness is, 
ae°a 2 . Equating this to surface energy (7a) and noting 
that the capillary pressure is linearly related to the stress, 
gives the critical capillary pressure for opening the crack, 

1 /3 

~ A (f ) 2/3 (^j , where '7' is the surface ten- 
sion of the solvent, R is the radius of the particles, and A 
is a constant. The objective of this paper is to rigorously 
determine the value of A and investigate the consequence 
of this result. 



II. NUCLEATION OF CRACK 

Consider a colloidal packing saturated with water that 
is confined by solid boundaries at x± — ±L with free 
surfaces at x-i — x% = ±h. As water evaporates, the 
capillary pressure puts the packing in tension in the 
x\ direction so that it is free to contract along 22 and 
X3 (Figure 2). In this case, the strain is given by, 
£y = — e° (#j2^2j + ^i3<^3i)) where denotes the pre-crack 
values. Volume conservation over a unit volume of the 
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bed relates the strain to the particle volume fraction, 
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where, <p° is the volume fraction in the pre-crack state. 
The strain and stress fields are sought for a crack with 
extent —a<X2 < a in the plain stress formulation, 
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with the perturbed variables represented by the primed 
quantities and a << L, h. Substituting these in the con- 
stitutive equation (1) and retaining terms linear in the 
perturbed quantities gives, 
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where a bar over a variable implies a dimensionless quan- 
tity with stress and pressure rendered dimensionless with 
E = GM 0rc P £ _ rpj ie dimensionless stress for the pre-crack 
state is, 

= <*« [-P° - 2(e°)] - 4(e°) (fe<% + <y a «3j) • (5) 

Since we consider the plane stress case (ct 3 j = 0) and 
the bed is stressed only in the X\ direction, P° — — 6e°, 
<j°h = 4e° and <j° 2 2 = 0- 

The total amount of particle phase remains constant in 
the packing, and so the particle volume fractions before 
and after cracking are related, 
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two cases, cr'ij — C^ef* where efj is the engineering strain 
and C is the stiffness matrix. The latter is given by, 
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in the short-time limit and, 
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in the long-time limit. While the original constitutive re- 
lation in the pre-crack state (1) is for an isotropic solid, 
(8) and (9) imply that the relaxation resulting from the 
presence of the flaw under the imposed conditions is that 
for an anisotropic solid. A similar observation is noted 
by Russel and co-workers for the more accurate consti- 
tutive relation based on the Hertzian contact mechan- 
ics. For convenience, we write the constitutive equa- 
tions as, Ai = SijTij where both A and S are 6-by-l 
column vectors. The components of A and S are given 

b y> e n) e 22) e 33) e 23> e 3i! e i2 and a ' n , (^'22, ■-, c'i2 respec- 
tively, and the 6-by-6 coefficient matrix, S, is the com- 
pliance matrix. The elements of S are determined from 
C. Note that the present case relates to the case of a or- 
thotropic anisotropy with plane stress condition. There- 
fore, S has only seven non-zero elements. The displace- 
ment field in case of plane strain is easily obtained using 
the procedure outlined in the next section except that 
the components of the compliance matrix for the plane 
stress problem (SV, ) are replaced with, 

n _ q _ S i3 S j3 
"-533 

for the plane strain case. 



(i,j = l,2,...,6) 



III. ASYMPTOTIC ANALYSIS NEAR A 
CRACK TIP 



where, 4> rcp is the random close packing and the strain 
is taken to be zero when (f> = 4> rcp . The time evolution 
of stress and strain around the crack can be further sub- 
divided into two limiting cases [3] i.e. the short time 
and the long time limits. In the short time limit, the 
impact of crack formation on the stress and strain varia- 
tion is such that it would occur in the absence of solvent 
flow, suggesting that the material will be incompressible. 
Thus, in the short time limit and for e° <C 1, (6) reduces 
to, 
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Since we shall consider only the plane stress case here 
(ct'3 3 =0), P' = e 22 + 7e 33 . At longer time scales, liquid 
flows so as to eliminate pressure variations, giving us the 
required condition for the long time limit, P' = 0. The 
perturbed stress and strain are compactly related in the 



The knowledge of the stress fields in the neighborhood 
of the crack tip is essential in determining the strength 
of the packed bed. Since the perturbed stress is linear in 
perturbed strain, we draw upon the mathematical tech- 
niques developed in the solid mechanics literature to de- 
termine the stress field near the tip of a crack [21-25]. 
The coordinate system (xx,x 2 ) for this analysis is shown 
in Figure 2 where the origin is placed at the crack tip so 
that x\ = Xi, x 2 — X2 + a and the variables have been 
rendered dimensionless using the characteristic length of 
the solution domain. The momentum balance equation 
in x\ and £2 directions in the absence of body forces are 
given by, 
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Following Sih et al.[21] and Hoenig[26], we related the 
stresses to stress correlation functions, x and ip, via 
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with eij and 5ij being the second order alternating and 
Dirac delta tensor respectively, and i,j allowed values of 
1 or 2. Note that (11) automatically satisfies (10). The 
above expressions along with the constitutive relation are 
substituted in the compatibility equations, 
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8x28x3 dx^dxi 8x 2 8x18x2 
to give, respectively, 

5'iiX,2222 + (25*12 + 5 l 66)x,ii22 + 522X,im = 0, and 

(13) 

5 4 4</\ll + ^55^,22 = 0. (14) 

The degree of anisotropy in the material can be judged 
by rewriting (13) differently, 



V 4 x + <5ix,im + <5 2 x,2222 = 
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where 1+Si = Q 2 f % 2 Q and 1 + So = s 2 f"e are indica- 



tors of anisotropy in the material. The difference in the 
values of 5i and 82 originate from the fact that the bed 
is held along '1' direction and perturbed along '2' and 
'3' directions leading to a directional perturbation of the 
stress field. When \5i\ -C 1, X satisfies the biharmonic 
equation, i.e. the material is isotropic. In the current 
problem, Si are —0.12 and 0.33 in the short-time limit 
and —0.24 and 1.4 in the long-time limit suggesting that 
the anisotropy is significant and cannot be ignored. 

(13) and (14) are a pair of decoupled equations in x 
and ip, 

U X = (16) 
L 2 tp = (17) 

where the differential operators are given by, L2 = 

2(5*12 + ^66) - 
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522 ^4 ■ Lekhnitskii[22] has shown that L 2 and L4 can be 

decomposed into two and four linear operators of first or- 
der respectively of the form D k = 8/8x2 — ^8/8x1 such 
that DiD 2 D 3 D4X = and D 5 D e tp — 0. Substitution of 
D k in L4 and L 2 shows that fj,k are roots of the polyno- 
mial operators, I4 = Sunt + (25i2 + 5 6 6)a4 + S22 = 
and I2 = S^fi 2 . + S44 — 0. Then, the stress correlation 
functions can be written as, 
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where the bar on fj,i represents the conjugate complex 
number. Further, Lekhnitskii[22] has shown that for the 
elastic energy of the packing to be positive, the roots 
cannot be real. Consequently, the general expression for 
the stress function will involve the real part of both the 
complex conjugates, 
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where Zi = Xi + HiX2 and Re represents the real part. 
Since the stresses are related to the derivatives of the 
above functions, it is convenient to assume the following 
functional form, 
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so that the stresses are given by, 



= G k {z k ) for ke [1,2] 



(20) 



o-'h 



c 22 



C 12 



O 31 = 



f'23 



d 2 XR 
8x\ 

d 2 XR 
8x\ 

d 2 XR 
8x18x2 

8x 2 
dtpR 



= 2 Re 



^ 2 dGi 
1 dzi 
dGi 

dzi 



= 2 Re 

= -2 Re 
= 2 Re 



Hi 



dG, 



dxi 



M3 



2 Re 



dz 3 
~dG 3 



M2 

dG 2 ' 

dz 2 _ 
dGi 
d^i 

, and 



,dG 2 
dz 2 



M2 



dG 2 

dz 2 



dz 3 



(21) 



We can also determine the strain components in terms of 
the stress function, for example, 
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Integrating the above expression with respect to Zj gives 
the displacement field in the '1' direction, 
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where, py = Sufi 2 + Si 2 - Si 6 (J.j- 

Following a similar procedure for the remaining strain 
components, all displacements are determined, 
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FIG. 3. Shifted coordinate system for the asymptotic analysis 



where 



displacement components near the crack tip become, 
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Note that the perturbed normal stress at the crack sur- 
face for the current problem has a finite non-zero value 
(— er n) which will require minor modifications to some 
of the above expressions and is dealt with towards the 
end of this section. The stress intensity factor (K) for 
mode I crack is defined as, 
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Similarly, K 2 and K 3 can be obtained from cr' 12 and a'13. 
Thus the stress intensity factor for the three modes can 
be written compactly, 
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Next, we ascertain the functional form of Gi. Rice [23] 
has shown that the J-integral, 

J = [ ( Wdxi — a'ijrij-^AS ] , has the same value 

J dp V dx 2 J 

for all integration paths surrounding crack tips in two 
dimensional fields of linear or nonlinear elastic materials. 
Here, W is the strain energy density, n.j is the normal 
to the chosen path, and S is the distance along the path 
dfl. Assuming ^j 1 oc zf in the neighbourhood of the 
crack opening, W oc z 2v and c'ij^jgf^ ~ z i 2p ■ Hence, 
J ~ Zi 2p+1 . Since the value of J should be independent 
of the path, p — — i. Thus, we assume G = Biyjlazijix 
for a flaw of size a where Bi is the stress function am- 
plitude. For a stress free crack surface, the stress and 
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The perturbed displacements of the crack surface can be 
found in terms of the distance from tip along the crack 
surface, 9 — > |, r = ( 
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where Q" 1 = Im{pN _1 I~ 2 }. 

The above analysis gives the functional form of the 
stress and strain fields close to the crack tip in terms 
of the unknown K. In order to determine the stress in- 
tensity factor, we assume a finite sized crack with an 
elliptical shape such that the minor axis of dimension- 
less length 2c is small compared to the major axis (2a), 
a = r « 1. Eshelby [27] has shown that for an elliptical 
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crack, the strain is uniform around the crack. Following 
Hoenig[26], the displacement of the crack surface is given 
by ($ = 1,2,3), 



U'i = A U X! = # 
and the strains by, 
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Thus, the perturbed stress at the crack face for this sim- 
ple crack model is related to a'if. — Cufa- Here, the 
origin of the coordinate system (xi,x 2 ) lies at the center 
of the ellipse with x 2 directed along the major axis (Fig- 
ure 2). Writing the crack face displacement in terms of 
the coordinate system with the origin placed at the crack 
tip (Figure 3), 
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where £ is the distance from the crack tip along x 2 . Com- 
paring (30) with (31) we get 
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which relates the unknowns, K and [0]. In order to com- 
plete the problem, we determine the elastic energy re- 
leased from the simple crack model, 



and equate gl = 2 J giving 
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where J is value of the standard J-integral [23], deter- 
mined in the limit as the integration path is shrunk so as 
to lie along the crack face, 
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Comparing (32) and (34) gives the expression for in 
terms of the far field stresses, 

fa = Qr fc W- 

We are now in a position to write the elastic energy re- 
covered (dimensional) due to the formation of a finite 
length mode-I crack of length a, 



(36) 



The present problem requires the crack surface to have 
a normal stress ( —a°n ) and the far field perturbed 
stress to be zero. Consequently, the complex stress 
function G k {z k ) in (21) is replaced with Gk(zk) + T^Zfe 



(k = 1,2), where T k are real constants [21]. Substitut- 
ing the new expression for G k and applying the traction 
boundary condition at the crack surface, we get, 
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The above expression along with the definition for the 
stresses (21) suggests that both cr'n and cr'22 are influ- 
enced by the traction condition. However, neither o'yz 
nor along the crack face are effected implying that the 
energy calculations and the corresponding values of the 
stress intensity factors remain unchanged. 

The total energy for the system is £ = £+r, where T = 
4r/a corresponds to surface energy and 7 is the surface 
tension. Following Griffith's argument [25], the crack will 
be in equilibrium when, 
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which relates the far field stress to the crack length and 
surface tension. Since E is a linear function of the pre- 
crack strain which in turn is related to the far field stress, 
we have 
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A more useful relation is in terms of the capillary pres- 
sure, 



P°R\ /a\ 2 /3 



1/3 1 ^ GMcf> rcp R 



1/3 



(40) 

Thus the dimensionless critical capillary pressure is re- 
lated to the dimensionless crack length, 



(-P°)(a 2 / 3 ) = AW 1 ' 3 



(41) 



where, W = GM(f> rcp R/2 r y represents the balance of the 
elastic and surface energy and A is equal to 0.45 and 0.35 
for the short and long time limit, respectively. 



IV. NUMERICAL SOLUTION 

The stress and the displacement fields obtained in the 
previous section is applicable to regions close to the crack 
tip. For the full solution, the momentum balance equa- 
tions (10) for the perturbed stresses were solved numeri- 
cally for the control volume highlighted in Figure 2 using 
finite element method (DIFFPACK®). The boundary 
conditions are as follows, 

er'n = — a°i for x\ = 0, > x-2 > —a, 
u[ = for X\ = 0, —a > x 2 > — 1, 
u 2 = for < xi < 1, x 2 = 0, 
cr'22 = for < x\ < 1, x% = — 1, and 
er'n = for x\ = 1, > x% > —1 
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FIG. 4. (b) The basic grid before deformation, (b) The scaled 
surface deformation of the control volume for the short time 
limit. All displacements have been scaled with one-tenth the 
maximum displacement (u[(0, 0)). 



where a << 1 so that the stress and the strain fields close 
to the cracks are not influenced by the size of the control 
volume. 

Rectangular elements were used with adaptive refine- 
ment of the grid near the crack tip. The total number of 
nodes in the control volume were about 50,000. 
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FIG. 5. Variation of perturbed stresses for nondimensional 
half crack length, a — 0.3 and e° = 7.8 x 10 -4 : (a) Contour 
plot of er'n, (b) Gray scale plot of er'n close to crack tip, (c) 
contour plot of a'22, and (c) Gray scale plot of er'n close to 
crack tip. 



V. RESULTS AND DISCUSSION 

A. Numerical Simulation 

The perturbed stress, strain and displacements ob- 
tained from the numerical simulations are presented in 
Figure 4-6. Unless specified, all results pertain to the 
short time limit. Figure 4(a) and (b) presents the initial 
and deformed grid, respectively. The surface displace- 
ments in Figure 4(b) have been scaled so as to highlight 
their magnitude. As expected, the displacement at the 
center of the crack is maximum with the tip of the crack 
moving upwards. Since the perturbed stresses are zero 
at the control volume boundaries, the effect of the sur- 
face displacement at the crack faces can be observed at 
the boundaries. The equilibrium shape of the crack sur- 
face is elliptical and the ratio of the length of the mi- 
nor to major axis is very small (~ 10~ 3 ), both of which 
are in agreement with the asymptotic solution. Fig- 
ure 5 presents the simulated values of er'n and a'22 for 
the control volume. The contour plot (Figure 5(a) and 
(c)) and the gray scale plots of the region close to the 
crack tip (Figure 5(b) and (d)) demonstrate the sharp de- 
crease in stress with increasing distance from the crack 
tip. The contours are perpendicular to the symmetry 
surfaces (52 = and —0.3 < x± < 1,X2 = 0) as expected 
from the boundary conditions while they decay to zero 
at x~x = 1 and x 2 = — 1- Figure 6(a) and (b) present the 
perturbed pressure for the short time limit. Interestingly, 
the pressure is negative close to the tip suggesting that 
the solvent will flow towards the tip once the crack nu- 
cleates. This is borne out in the simulations for the long 
time limit (Figure 6(c) and (d)) where the particle con- 
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FIG. 6. Variation of perturbed pressure in the short time limit 
for a = 0.3 : (a) Contour plot of P' , (b) Gray scale plot of P 1 
close to crack tip. Variation of the particle concentration in 
the long time limit: (c) Contour plot of <f>, and (d) Gray scale 
plot of (f> close to crack tip. 



centration has reduced at the tip. Note that while the 
particle concentration should always be equal or more 
than the close pack concentration at all times, values of 
4> < fircp near the crack tip in the long time limit are not 
physical since no such constraint has been imposed in 
the simulation. Instead, extra solvent could accumulate 
at the crack tip between the crack faces. 
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FIG. 7. Displacement u[ for a — 0.3 along the crack face. 
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FIG. 8. Perturbed stress a'u along xi (x2 = —a) for a — 0.3. 



B. Comparison with Asymptotic Solution 

Figure 7 compares the spatial variation of perturbed 
displacement, u'x, along the crack face obtained from the 
simulation with that predicted by the asymptotic solu- 
tion. At the crack tip, = while for 1CP 4 < x\ < 
10~ 2 , the displacement varies as the square root of the 
distance from the tip. The disagreement close to and far- 
away from the crack tip is attributed to the limitation 
on grid refinement in case of numerical solution close to 
the tip (which is unable to capture the large variations in 
the stress) and to the non-applicability of the asymptotic 
solution far away from the crack tip. Figure 8 presents 

the spatial variation of a'u along x% away from the crack 

_ i 

tip for a — 0.3. As expected the stress diverges as x 1 2 
close to the crack tip and the prediction matches well 
with the numerical solution for 10 -4 < x\ < 10~ 2 . Far 
away from the crack tip, the perturbed stresses vanish. 
The angular distribution of stresses obtained from the 
asymptotic solution agrees with that from the numerical 
solution in Figure 9 at |f — <x| = 0.012. The distribution is 
somewhat similar to that obtained for the isotropic cases 
[28]. Figure 10 compares the angular distribution of the 
stresses at various radial distances from the crack tip, 
both in the short and the long time limits. The magni- 
tude of the stresses at a given location in the long time 
limit are lower than those in the short time limit. This 
decrease may be attributed to the flow of the solvent that 
relieves any pressure variation that develops in the short 
time limit. Compared to the short time limit, the angular 
variation of the stress in the long time limit show larger 
deviations from the isotropic case as also suggested by 
the values of Si in (15). 

The Griffith's criteria (41) shows that the pressure re- 
quired to open a crack increases with decreasing size 
of the crack. Since the maximum dimensionless cap- 
illary pressure is about [29] 5.3, the largest allowable 
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FIG. 9. Angular variation of perturbed stresses for short time 
at |f-a| = 0.012 for a = 0.1. 



flaw which will not crack the sample in the short time 
limit is, a max = 0.025\/VF which suggests that packings 
containing particles of larger size and/or higher shear 
moduli can resist cracking more effectively. Recently, 
Tirumkudulu and Russel[3] have derived the expression 
for the critical capillary stress to drive an infinite crack 
through a drying colloidal thin film bound to a substrate, 
(-P^)(h 2 / 3 ) = 0.23W 1 / 3 . Comparing the critical capil- 
lary pressure for the two cases suggests that when a <C h, 
a significantly larger capillary pressure is required to ex- 
pand a finite flaw in the film compared to that required 
to drive an infinite crack, 




These results are in line with the recent theoreti- 
cal results obtained by Russel et al.[18] using the more 
accurate constitutive relation and experiments measur- 
ing the critical capillary pressure for various particle 
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FIG. 10. Angular variation of perturbed stresses at different 
radii (R = f — a), (a)-(c) in the short time limit, and (d)-(f) 
in the long time limit for a = 0.3. 



VI. CONCLUSIONS 

We present the asymptotic analysis of the deformation 
field near a crack tip for a mode I crack in a two di- 
mensional colloidal packing saturated with solvent. The 
stress and strain fields are linearized about the pre-crack 
state to yield the stress intensity factor for the two di- 
mensional elastic field which is then related to the surface 
energy using the well known Griffith's criterion for equi- 
librium cracks. The calculated quantities are then com- 
pared with the numerical solution for the full problem. 
The main findings can be summarized as follows: 

• Perturbation in the displacement and stress field 
due to the presence of crack introduces anisotropy 
in the material which can be quantified by (15). 

• The stress and displacement fields close to crack 
tip are given by (26) where the expression of B is 
obtained from the components of stiffness matrix, 
C. 

• The critical pressure required to open a flaw of 
length 2a varies inversely with the crack length to 
the two thirds' power, 



-P° = A(GM^ rcp ) 1/3 



2/3 



packings [19]. 

The energy release rate, Q = 2 J, is related to the stress 
intensity factor through the standard relation, Q = -j^ 
where the effective elastic modulus for the packing, 

E eS = -^, (42) 

accounts for the particle size and packing, and also for 
the anisotropy resulting from the nucleation of crack. 

It is important to note the limitations of the analysis 
presented here. The boundary condition a'n = —<j°i at 
the surface of the crack implies perturbed strains of the 
order e° close to the crack surface, which is inconsistent 
with the linearization in (4). The same applies to the 
diverging perturbed strains at the crack tip as predicted 
by the linear analysis. The extent of errors introduced by 
such approximations can be accurately determined only 
by solving numerically the full non-linear momentum bal- 
ance equations. However, the recent experimental evi- 
dence of diverging stresses close to the tip of a crack in a 
colloidal packing supports the overall trend predicted by 
the linear analysis. 

Finally, the analysis presented here is general, in that 
the results relating to the asymptotic forms of the stress 
and displacement component, and the related expression 
for the energy release rate can easily be obtained for any 
other constitutive equation for a saturated packed bed 
once the stiffness matrix for the linearized equation is 
known. 



where A is equal to 0.45 and 0.35 for the short and 
long time limit, respectively. It is independent of 
the particle size. 

• The maximum flaw size that can resist cracking and 
result in a crack-free packing is set by the maximum 
possible capillary pressure, 



.4 
5^3 



3/2 / GM(f> rn „R 3 



u rC pi 

17 



1/2 



Colloidal beds containing large particles with high 
shear modulus are less susceptible to cracking. 

• When a <C h, the critical capillary pressure re- 
quired to expand a flaw is much larger that that 
required to drive an infinite crack in a film of thick- 
ness, h. 
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